Theory of light diffusion in disordered media with linear absorption or gain 
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We present a detailed, microscopic transport theory for light in strongly scattering disordered sys- 
tems whose constituent materials exhibit linear absorption or gain. Starting from Maxwell's equa- 
tions, we derive general expressions for transport quantities such as energy transport velocity, trans- 
port mean free path, diffusion coefficient, and absorption/gain length. The approach is based on 
a fully vectorial treatment of the generalized kinetic equation and utilizes an exact Ward identity 
(WI). While for loss- and gainless media the WI reflects local energy conservation, the effects of 
absorption or coherent gain are implemented exactly by novel, additional terms in the WI. As a 
result of resonant (Mie) scattering from the individual scatterers, all transport quantities acquire 
strong, frequency-dependent renormalizations, which are, in addition, characteristically modified by 
absorption or gain. We illustrate the influence of various experimentally accessible parameters on 
these quanitities for dilute systems. The transport theory presented here may set the stage for a 
theory of Random Lasing in three-dimensional disordered media. 
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I. INTRODUCTION 

Despite its long and venerable history, light propaga- 
tion in disordered media continues to be a fascinating and 
intensely studied topic. In particular, the discovery of the 
coherent back scattering peaki 2 ^ has triggered an inten- 
sive search for Anderson localization 4 of light. Clearly, 
the unambiguous demonstration of Anderson localization 
of electromagnetic radiation in an appropriate disordered 
dielectric medium in which dephasing and interaction ef- 
fects can be neglected, would mark a major triumph of 
our current understanding of wave propagation. 

In fact, early works reported anomalously low val- 
ues of the diffusion constant in strongly scattering 
suspensions^. However, it was soon realized^**! that 
these low values of the diffusion constant are associated 
with the resonant (Mie) scattering of individual scat- 
terers. This leads to a frequency-dependent dwell-time 
that has to be added to the time-of-night between suc- 
cesive scattering events. As a result, the energy trans- 
port velocity acquires a corresponding renormalization, 
while the transport mean free path remains essentially 
unchanged^. Similarly, an analysis of the coherent 
back scattering peak together with the dependence of the 
transmission through disordered semiconductor powders 
on the sample size, suggested a scaling behavior«2> consis- 
tent with the onset of Anderson localization^. However, 
a reexamination of these data ignited a heated debate as 
to how to discriminate Anderson localization from ab- 
sorption cffcct q 11 ! 12 . Subsequent experiments on sim- 
ilarly strongly scattering semiconductor powders 1 ^ did 
not produce evidence of Anderson localization and the 
experimental data could be well explained using a re- 
cently developed effective medium theory that incorpo- 



rates the resonant scattering effects alluded to abov o 14 ' 15 . 
This unsatisfactory state of affairs has generated renewed 
interest in determining novel and unambiguous pathways 
to Anderson localization of light. 

One class of highly interesting systems for multiple 
scattering of light are disordered dielectric media whose 
constituent materials exhibit one or more forms of op- 
tical anisotropics. Moreover, most optical anisotropics 
exhibit a certain degree of tunability through external 
control parameters, thus creating the possibility of a tun- 
able disorder. For instance, Faraday-activty in multi- 
ple scattering systems breaks the time reversal symme- 
try of Maxwell's equations, leading to profound modifica- 
tions of the coherent back scattering peakiSiiLi^ and in 
transport the optical analogue of the Hall-effect has been 
observed 1 ^. Similarly, disordered nematic liquid crys- 
tals exhibit anisotropic coherent back scattering^ and 
anisotropic light diffusio n 21 i 22 . However, the strength of 
multiple scattering in bulk nematic liquid crystals is de- 
termined by the fluctuations of the nematic director field 
and the difference between the liquid crystal's ordinary 
and extraordinary index of refraction. For an observa- 
tion of Anderson localization of light in such systems, 
it will become necessary to enhance the multiple scatter- 
ing effects, for instance, by infiltrating the nematic liquid 
crystal into the void regions of strongly scattering pho- 
tonic crystals^. 

Equally intriguing is to combine the effects of multi- 
ple scattering of light with optical gain and to investi- 
gate how the two phenomena mutually influence each 
other. From the optical gain point of view, diffusing 
light will spend more time interacting with active ma- 
terial in a characteristic volume than ballistically propa- 
gating light. If this interaction time exceeds the (sponta- 
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neous) decay time of the active material, avalanche-like 
intensity bursts, induced by incoherent feedback, could 
occur. In fact, early theoretical work2i suggested this 
very possibility and has recently been observed^. From 
a wave propagation point of view, optical gain increases 
the relative weight of long trajectories in the sample and, 
therefore, will modify a number of wave interference ef- 
fects such as coherent back scattering 26 and possibly An- 
derson localization. Very recent experiments point to the 
possibility that such long trajectories provide a feedback 
mechanism which leads to modes with narrow laser-like 
emission lines that extend across the entire sample 2 ^. 
However, the coherence properties of the emitted light 
still need to be analyzed in order to establish that laser 
action is indeed taking place in these systems. True laser 
action from localized regions in disordered dielectric sam- 
ples with optical gain, termed Random Lasing, has been 
observed2Si22*^ earlier, where measurements of the pho- 
ton statistics of the emitted lighted have unambiguously 
demonstrated a coherent laser feedback mechanism. It 
has been discussed in the form of varying degrees of cou- 
pling between so-called quasi-states^S. For a recent re- 
view of multiple scattering in amplifying media, we refer 
the reader to Refi^i. Random Lasing has potential ap- 
plications ranging from micro-lasers and optical finger- 
print markers^ to the detection of cancerous tissue^ 6 .. 
In addition, we want to note that recently, electrically 
pumped Random Laser action has been achieved in Nd- 
doped powders2L2&, thus creating potential applications 
of these systems in omnidirectional lighting devices and 
displays. 

Despite these exciting developments, a convincing con- 
nection between Random Lasing and Anderson localiza- 
tion of light has not been demonstrated as of yet. This 
may be attributed to the fact, that to date the the- 
ory of random amplifying media either employs purely 
numerical methods in one spatial dimension 3 — or treats 
the multiple scattering part within certain approxima- 
tion schemes that cannot account for the interference ef- 
fects that lead to Anderson localization. These schemes 
include modeling the electromagnetic wave propaga- 
tion through diffusion equations for the intensity^Mii or 
within the so-called ladder-approximation of the Bethe- 
Salpeter equation^. As a result, it is unclear, whether 
Anderson localization of light is a necessary condition for 
true Random Lasing (coherent feedback) nor whether the 
modified coherence properties in the lasing state have, in 
turn, an influence on the transition from the diffusion 
regime to the Anderson localized regime itself. In fact, 
the very concept of describing the Anderson localization 
transition in terms of a vanishing diffusion coefficient as 
an order parameter, familiar from systems with energy 
or particle number conservation, becomes questionable 
in dissipative or active media, where another channel for 
change of the energy density exists besides diffusion. 

In the present paper, we report our progress towards 
answering these questions. We develop a fully vectorial 
transport theory for multiple scattering of light in ran- 



dom media whose constituent materials are isotropic and 
exhibit linear absorption or gain. In Section II, we de- 
rive the tensorial kinetic equation for the intensity cor- 
relation function of electromagnetic radiation. The con- 
servation of energy in media without loss or gain is in- 
corporated in a field theoretical way by means of an ex- 
act Ward identity (WI), the effects of loss, gain as well 
as frequency dependence of the material parameters be- 
ing represented by additional terms. The proof of this 
generalized WI is presented in Section III. Subsequently, 
we solve the kinetic equation in the hydrodynamic limit 
by formulating, with the help of the WI, the continu- 
ity equation in Section IV and Fick's law in Section V, 
respectively. These equations relate the energy density 
correlation tensor and the energy current correlation ten- 
sor to each other. They allow to identify, in the hy- 
drodynamic limit, exact expressions for quantities like 
the energy transport velocity, the transport mean free 
path and the absorption/gain length in terms of the ir- 
reducible single-photon self-energy and the two-photon 
irreducible vertex. We recover the well-known Mie scat- 
tering renormalization term to the transport velocity^, 
albeit in vectorial form, along with additional, charac- 
teristic renormalizations originating from absorption or 
gain. Section VI features numerical results for the en- 
ergy velocity in dilute systems of spherical scatterers, 
where all quantities can be evaluated within the inde- 
pendent scattering approximation. The Mie resonance 
dips in the energy transport velocity acquire a charac- 
teristic, absorption-induced broadening and gain-induced 
narrowing, which may be interpreted as due to the reduc- 
tion (enhancement) of higher-order scattering contribu- 
tions by absorption (gain). In contrast to previous ap- 
proaches, our theory can be extended to systematically 
include wave interference effects ("Cooperon" contribu- 
tions) and, thus, to address the Anderson localization 
transition by virtue of a self-consistent extensior&M44i2£ 
together with a replacement of the linear absorption/gain 
through a coupling to the master equation of the active 
medium^. 



II. BASIC THEORY OF ELECTROMAGNETIC 
WAVE TRANSPORT 

We consider the propagation of light in a system of 
randomly positioned scatterers with isotropic dielectric 
constant e s immersed in a host medium with isotropic 
dielectric constant eh- In addition, we allow for the pos- 
sibility of having absorption or amplification in both the 
scatterer material and the host medium by ascribing an 
imaginary part to the respective dielectric constants. By 
virtue of the Kramers-Kronig relations between real and 
imaginary parts of the dielectric constant, we are required 
to consider materials with complex frequency dependent 
dielectric constant e s = e s (uj) and eh = £h(u). The re- 
sulting dielectric constant 

e(f,u) = e h + (e s -e h )V(r), (1) 
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describes the arrangement of scatterers through the func- 
tion V(r) = ^2^Ss{r — R) which consists of a set of 
localized shape functions Sg(r) of the invididual scat- 
terers at random locations R. The corresponding wave 
equation for a time harmonic wave with frequency lu and 
amplitude E^if) reads 



V x V x E u {f) - -^-e(r,w)E u (r) =uj J u {f\. 



(2) 
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FIG. 1: Diagrammatic representation of the Bethe - Salpeter 
Equation, Eq. indicating the tensorial structure as well 
as the full index notation. 



Here, c denotes the vacuum speed of light, and J w (r) = 
(4iri/c 2 ) j u (f) represents an excitation of the wave field 
through an external current source j. The dielectric 
function e(r, lu) carries a positive or negative imaginary 
part for absorbing or coherently amplifying media, re- 
spectively. 

For a specific realization e(r, lu) of disorder, the formal 
solution to Eq. is given in terms of the Green's tensor 
(of 2nd rank in the spatial vector components) G(r, f, lu) 
as 



E^f) = / d 3 r'G(ry,uj),U7<). 



(3) 



For the analytical developments as well as for the com- 
parison with experiments it is advantageous to calculate 
disorder-averaged quantities. Since these are translation- 
ally invariant, the transport theory must be formulated 
for correlation functions rather than single-particle prop- 
erties. Denoting the average over disorder realizations 
by (...), and changing to a Fourier representation of 
averaged quantities, the disorder averaged Green's ten- 
sor G u (r — r 1 ) = (G(r, r 1 ,lu)) and its Fourier transform 
Gg may be expressed in terms of the free Green's tensor 

Go(k, lu) in the background medium, 

G (k, lu) = (e h {Lu){u 2 /c 2 ) - |fc| 2 ) ~* (l - kk T ) , (4) 

and the self-energy tensor Eg, which represents the ef- 
fects of scattering by the random perturbation or "scat- 
tering potential", {lu 2 /c 2 )(e s (w) — €h(u)))V(r), as 



(5) 



In the above expressions, we have introduced the three- 
dimensional unit tensor 1 and the dyadic product kk 
of unit vectors k in the direction of k. Throughout 
this paper the propagators G are understood as the re- 
tarded ones and complex conjugated propagators G* as 
the advanced ones. For practical calculations, the self- 
energy tensor Eg has to be evaluated within consistent 
approximations such as the independent scatterer ap- 
proximation (see section VI) or the Coherent Potential 
Approximation^^. 

The disorder-averaged field correlation tensor is 



dcfinied as 
I{r u r 2] L0 X ,LU 2 ) = (E ul (n)E^(r 2 )}) 

= J d 3 r 3 J d 3 r A (fj., f 2 , r 3 , r 4 ) 



xS Ul , U2 (r 3 ,f 4 ), 



(6) 



where the system's (averaged) correlation tensor 
(of fourth rank) T ull!UJ2 (fi, f 2 , r 3 , f 4 ) is independent 
of the source correlation tensor S u!l _ u!2 (r 3 ,r 4 ) = 
Juji {^i)J^ {^i) an d is given in terms of the disorder av- 
eraged tensor product of (unaveraged) Green's tensors 
according to 

T Ul ,u 2 {ri,r 2 ,rz,ri) = 

(G(f 1 ,f 3) a;i)«)G*(f 2 ,f4,W2)) , (7) 

where (. ® .) denotes the tensor product of two 2nd rank 
tensors operating in the space of retarded and advanced 
propagators, respectively. Similar to the disorder aver- 
aged Green's tensor G|, we introduce the spatial Fourier 

transform r^=» (Q, fi) of the correlation tensor 

r^-, (Q,fi) = (G u+ (?+,? + ) g:_ (?_,<?_)>, (8) 

where the transition to center of mass and relative fre- 
quencies, fi, lu, and momenta, Q, q, respectively, with 
^1,2 = lu ± fi/2 = lu±, qi t2 = q ± Q/2 = q± and 
Q3.4 = <f i Q/2 = q±, facilitates an investigation of the 
correlation tensor's long-time (fi — ► 0) and long-distance 
(\q\ — ► 0) behavior, fi, Q are associated with the time 
and position dependence of the electromagnetic energy 
density in the system, while lu represents the frequency 
of light . 

(Q, fi) can be expressed in terms of the irreducible 

vertex tensor 7—, (Q, fi), the two-photon analogue of the 
self-energy tensor, via the Bethe-Salpeter equation 



r&,(Q,n) = G£ 



[(27r) 3 S(q - $01' 



d 3 q" 



7 ^(Q,fi)r^ v (Q,fi) 



(9) 



In Eq. both the irreducible vertex tensor 7~»(Q, fi) 

and the correlation tensor (Q, fi) are tensors of fourth 
rank which operate both in retarded and advanced space 
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(seeEq. ©). The notation 7^ (Q, ft) (Q, ft) in Eq. 

implies contraction of both, retarded and advanced, 
intermediate indices, 



7&/(Q,fi)r£ 4 .(Q,n) 



a/3 err 



(10) 



where summation over repeated indices is implied, com- 
pare Fig^ For tensor products of 2nd rank tensors in re- 
tarded and advanced space, B,C and E,F, respectively, 
we have the multiplication rule 



(B ®E) (C <S> F) — (BC) (8 (EF), 



(11) 



where BC and EF denotes the standard (matrix) prod- 
uct between tensors of second rank. It follows the general 
identity 



(g£ : )*=(gS+ 



-1 



1 



(<*)*) 



(12) 



which allows us, after integration over the momentum <f, 
to rewrite the Bethe-Salpeter equation © as a kinetic 
equation (a generalized Boltzmann equation) for the sys- 
tem's integrated intensity correlation tensor 



d 3 q' 



(2tt) 3 qq 

Explicitly, the generalized Boltzmann equation reads 
-Ag u l ® 1 - AL g -+ A5 

d 3 g 



(13) 
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AG£ 



10 1- 



(2tt 



3 Iqq" ^q" 



(14) 



Here, we have introduced short-hand notations for cer- 
tain differences and, for later use, sums of 4th rank ten- 
sors, 



Ag w (ft) 


= g(v+) 






(15) 


AL q iQ) 


= Hq+) 


® 1 - 1 c 




(16) 


A^(Q,fl) 


= « 


5 1-1® 


(*£)' 


(17) 


AG%(Q, fl) 




5 1 - 1 (g) 




(18) 


US%(&,fl) 


- m « 


5 1 + 1® 




(19) 


IlG%(Q;fl) 




5 1 + 1® 




(20) 



which employ the definitions g(oS) — (w 2 /c 2 )e/ l (w) and 
L(r) = V x Vx with corresponding Fourier representa- 
tion L(q) = —q x qx. 

The physical interpretation of the generalized Boltz- 
mann equation, Eq. I|14fl starts with the long-time limit, 



i.e. small fl expansion of the first term on the left-hand 



side (l.h.s.), Ag UJ (fl) = ig t 



(0) 



where 



7 (0) 



,(!) - 



2^Im [e h («)] 



Re 



9 w -^-eh.(w) 



(21) 
(22) 



It is seen that the term linear in fl, corresponds 
to a rate of change of the correlation tensor (Q, fl) 
(term of 0(fl)), while the term of O(ui ) describes ab- 
sorption or emission by the host medium. The latter is 
non- vanishing only if Ime^a;) ^ 0. Similarly, the second 
term on the l.h.s. represents a drift term (first order in 
q) for fl). The third term on the l.h.s of Eq. ((Hj) . 

AE^(Q, fl), embodies single-particle scattering from the 
external random perturbation, while the terms on the 
right-hand side (r.h.s.) represent an effective two-particle 
collision integral induced by disorder averaging. 

The physical interpretation of the generalized Boltz- 
mann equation suggests the subsequent strategy for ob- 
taining the solution in the hydrodynamic limit (fl — » 0, 
Q — > 0), where the collective mode dynamics are gov- 
erned by the conservation laws of the system. The elec- 
tromagnetic wave equation (J2J) being 2nd order in time 
implies (in a loss- and gainless medium) the local conser- 
vation of the energy density rather than the intensity^. 
Consequently, we will seek solutions for the energy den- 
sity correlation tensor P and the energy current density 
correlation tensor J, and recast the kinetic equation in 
terms of these quantities. This is possible because, as 
seen in section IV, P and J are the leading coefficients 
of the correlator fl) in a small Q, fl expansion. 

Before we develop this solution in Sections IV and V, 
we will derive in the next section an exact WI for vec- 
tor waves, which relates the photon selfenergy £ and the 
irreducible two-photon vertex gamma to each other. It 
embodies local energy conservation as well as disspation 
or gain induced deviations in a field theoretical way. 



III. WARD IDENTITY 

The derivation of the WI for disordered electronic sys- 
tems has been demonstrated for the first time by Voll- 
hard and W61fl o 43 i 44 using a diagrammatic approach. A 
corresponding WI for scalar classical waves has been de- 
rived, using an algebraic approach, by Barabanenkov et 
al4& and, using a diagrammatic technique, by Kroha et 
&\&. Barabanenkov et al. have generalized their deriva- 
tion to electromagnetic waves^S^i. Subsequently, the 
correct form of this WI has been the subject of heated 
debatei*2iSi£4, where a consensus has been reached in Ref. 
IHil However, to date, frequency dependent and/or com- 
plex dielectric functions have not been included in the 
derivation of WI for classical waves. Neither have all the 
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implications of the WI on the renormalization of trans- 
port quantities been discussed. 

Therefore, the aim of this section is to derive a WI 
for electromagnetic waves in the presence of frequency 
dependent, complex dielectric functions. In the following 
proof, we choose to follow the approach by Barabanenkov 
et aliSi. A detailed discussion as to how the WI affects 
the various transport properties is postponed to sections 
IV and V. 

We start from the Green's tensor before impurity av- 
eraging whose equation of motion (see Eqs. an d ©) 
we write as 

(i(ri) + <?M - h{u!)V(fij) (23) 

xG(fi,f 2 ,uji) = 5(n - 7=2)1. 

Here, we have introduced the, in general, complex quan- 
tity 



h(u) = -(e s (u)-e h (u))(uj 2 /c 2 ) . 



(24) 



Multiplying Eq. (|23|) with h(uj 2 )*G (^3,^4,^2) within 
the appropriate tensor subspace yields 

h*{u 2 )L{rx)®\ + g{(Ji)h*{u 2 )l®l (25) 

-h{ui)h*(LJ 2 )V{ri)t®l 

xG(fi,f 2 ;w 1 )<g>G (r 3l r 4 ;u2) 
= h*{u 2 )6(fi - r^)!® G (7=3, ri;^)- 

In the time reversed case, i.e. starting with the complex 
conjugated equation, an analogous procedure leads to 

h(u 1 )l®L*(f f a)+g(cj2)h((jj 1 )l®l (26) 

-h(u)i)h*(w 2 )V(?3)l®l 

xG{fi,f 2 ;vi)<g>G* (f=3, f 4 ; w 2 ) 
= h(uji)5{r 3 - fi)G{ri,f 2 ;w\)®l. 

Upon substracting Eqs. (|2*5|) and Eq. {25)1. followed by 
averging over disorder and taking the limit of r\ — > f 3 , 
we combine retarded and advanced quantities to what 
Barabenenkov et ali^i refer to as the pre- Ward identity. 







lim 

fi-»f 3 



fc*(w_).L(ri)<8>l - h{uj+)l®L (f 3 ) 

(27) 



- \g(w+)h*{u-) - g*{tjj-)h{u + )] 1®1 

(G(fi, r 2 ; uj+) <g> G (f 3 , r 4 ; 
-(^(w-)^^ -f 2 )l(g>G*(f 3 ,f 4 ;to-) 

-h(u+)5(f 3 -f 4 )G(fi,f 2 ;cJ + ) ®1) I 



The final step in obtaining the WI is to transfer the 
pre- Ward identity to momentum space and to cancel or 



simplify several terms with the help of both the Bethe- 
Salpeter equation and the generalized Boltzmann equa- 
tion (for details we refer to the work of Barabanenkov et 
al^i). 

This finally yields the Ward identity for electromag- 
netic waves with complex frequency dependent dielectric 
functions. 



AE£(Q;fi)- 



d 3 c 



; J (2tt) 3 
h{uij r ) + h*{ujJ) 

dy 



IIE£(Q;f>) (28) 



(2nf 



UG^(Q;n)^(Q,n) 



As compared to the case of electronic wave propagation 
in a disordered solid^iMi, the non-zero r.h.s. of Eq. i|28|) 
constitutes a novel term which originates from the de- 
pendence of the "scattering potential" (uj 2 /c 2 )(e s (uj) — 
eh(u>))V (r) on the light frequency as well as from its 
possible imaginary part, and will lead below to a renor- 
malization of the energy transport velocity^i^. In fact, 
in the case of frequency independent and real dielec- 
tric constant, we have that the prefactor of this term, 
(h(u) + ) ~ h*(uj-))/{h{uj + ) + h*(ui-)) = 4wf7/(4w 2 + ft 2 ), 
takes on a form that has been discovered earlier—. Since 
in this case the r.h.s. of Eq. I|28[) is oc f2 it leads to a renor- 
malization of the energy transport velocity 8 . However, in 
the present case of absorptive or amplifying media, the 
r.h.s. contributes also in zero-th order in Of2, signal- 
ing that absorption and gain induce more severe effects 
than renormalizing the energy transport velocity, namely 
a mass term for the diffusion modes, see below. 



IV. CONTINUITY EQUATION 

We now proceed with solving the kinetic equation i|14fl . 
From the 2nd order wave equation J2J it follows that for a 
homogeneous system (e(r) = tu — const.) the quantities 



Pe 



Je = i— #-(V£)-(V£+M 



(29) 
(30) 



obey the continuity equation dpg/dt + V • Je = 0, and 
may be interpreted as energy density and energy cur- 
rent density, respectively, where c — c/ yTiee^(w) is the 
phase velocity of the homogeneous medium. In order 
to obtain a similar relation for correlation functions in a 
random medium we, therefore, define the energy density- 
energy density correlation tensor P^, and the longitudi- 
nal energy current-energy density correlation tensor J" 



6 



in Fourier representation as 



p^(Q,n) = Lu 2 c- 



cl 3 



Q 



(2tt) 3 " 
d 3 



(31) 



J-(Q, n) = c^Hc- 1 !^ (g ■ Q) 0), (32) 



where 



(33) 



defines the averaged phase velocity tensor in the ran- 
dom medium^, and the energy transport velocity tensor 
v e (uj) is to be identified below. We stress that in Eqs. 
(131(1 and ((32(1 the products on the r.h.s. are contractions. 
For instance, in full index notation Eq. (|31|l reads 



a0 



■U? \c~ 2 ] 



(2nj 



3 Q 



ija/3 



where summation over repeated indices is implied. A 
corresponding expression can be written down for the 
current density tensor, Eq. (|32(l . 

The expressions Eqs. (|31|) . 132|) are the leading terms 
of a small fl, Q expansion of <6 45 i 47 , analogous to the 
so-called Pi-approximation to the standard Boltzmann 
equation 4 ^, 



AG£(0,0) [Pg(Q, n)A^ 
+(q-Q)J%(Q,n)A^ + 0(ft,Q 2 )(35) 
d 3 <jAG|(0,0) (36) 

d 3 q(q-Q) 2 AG%(0,0), (37) 



The tensor coefficients Ao and A\ above can be computed 
by projecting $^(Q,f2) onto its Oth and 1st moments 
with respect to the longitudinal current vertex (<f • Q), 
i.e. ontoPg(Q,fi) and J^(Q,iY) respectively. 

The continuity equation for Pe and Jg, can be derived 
by integrating the generalized Boltzmann equation, Eq. 
1(14(1 . over the momentum q and subsequent application 
of the WI, Eq. (|28|l . Upon considering the long-time, 
long-distance limit, we arrive at 



uj z c- z ld 3 qAG%{Q,tt) + -i fl (°)l® 1 + AH 



(38) 



The renormalization terms A(w) and A(w) originate from 
the nonzero l.h.s. of the WI, Eq. 1)28(1 . and are defined 
as 



A(w) = a (1) h+a (2) h 

(o) (i) 

AH = ^f^A^H, 



(39) 
(40) 



with the explicit forms of A^H an d A( 2 )(w) given by 
A (1) H = ~^/(0|nS|(O,O)AG^(O,O) (41) 

+ Jj0f UG ? (°. W*( > ^)ag|(o, o) j V 



3^(0,0) AGS (0,0) (42) 



dV 



3 [d n G£(0, 0)^(0, 0)AG£(0,0) 
f IIG^t (0, 0)[cVy^(0, 0)]AG| (0, 0) U^ 1 . 



(27T) 

(>) 3 



,(34) In Eqs. 1(38142(1 we have expanded the prefactor 



/i(w + ) + h*(uj-) 
to leading order in Q, where we introduced 



( 0) = Im [feH] 
w Rc [/i(cj)] 

(1) _ Rc [h(J)] Rc [d u h{u))] + Im [/i(w)] Im [d u h(u)] 
luJ ~ 2Rc [h(w)] 



(44) 



(45) 



Eq. I|38|) takes the form of a continuity equation by iden- 
tifying the energy transport velocity u B (w) such that the 
prefactors of the two terms on the l.h.s. are equal, and 
by multiplying with the inverse of that prefactor. In fact, 
this procedure defines the energy transport velocity, 

u E H = ^ y [c P r 1 Nl + AH]" 1 , (46) 

where the r.h.s. implies a contraction as explained in Eq. 
(|31|l . The continuity equation then reads, 



QJ% = 



(iy [1 ® 1 + AHP u 2 c; 2 d 3 qAG%(Q, n) (47) 

+ [18 1 + AHP 1 ^ 1 <g> 1 + A(w) 

We want to emphasize that in the long-time, long- 
distance limit, the continuity equation, Eq. 1(47(1 is ex- 
act. Also note that in the disordered medium transport 
occurs, with different velocities, both perpendicular and 
parallel to the polarization directions of the in- as well as 
of the outgoing wave field, since each scattering process 
alters the polarization. Accordingly, the averaged trans- 
port velocity v-e(lj) is a 4th rank tensor, as evidenced 
by Eq. 1(46(1 . Similarly, in a transport theory of a vector 
field the diffusion coefficient D(uS), the transport mean 
free path t T (u>), the absorption/gain length x a (u>), given 
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below by Eqs. I|55[l . Q56J1. and (|57|l . respectively, and other 
transport quantities are tensors as well. 

The microscopic derivation of the energy transport ve- 
locity, Eq. and the continuity equation, Eq. (|4"?)) . 
exhibit a number of important physical aspects: 

(1) The frequency dependent renormalization A (a;) of 
the energy transport velocity tensor consists of the "dwell 
time" renomalization A^(ui) alluded to in the introduc- 
tion which in the case of independent scatterers may be 
traced to their internal (Mie) resonances. 

(2) However, in the presence of a mismatch in absorp- 
tion or amplification ( "impedance mismatch" ) between 

the scatterer and the host medium (non- vanishing aL^), 
Eq. (|46|l features additional renormalization of the trans- 
port velocity. To the best of our knowledge, this renor- 
malization has not been discussed before. Numerical re- 
sults for the transport velocity in absorptive or active 
media will be presented in Section VI. 

(3) The r.h.s. of the continuity equation, Eq. (|47|l fea- 
tures as the first term the source contribution expected 
for correlation functions. We find that in the presence of 
either absorption or gain the second term on the r.h.s. 
does not vanish. Moreover, a "dwell time" renormaliza- 
tion A(uj) oc A' 1 ) (uj) occurs in the case of unequal absorp- 
tion or gain characteristics of host medium and scatterer 
(non- vanishing a^). While such a renormalization of 
either absorption or gain could be expected on physical 
grounds, our calculation constitutes the first quantita- 
tive analysis of this effect. Both absorption/gain related 
renormalization terms may have profound influence on 
lasing action in random dielectric media and further ex- 
plorations are necessary. 



FICK'S LAW AND DIFFUSION POLE 



In order to close the set of equations for P E and J E , 
another relation besides the continuity equation has to 
derived from the kinetic equation, Eq. I|14|) . which relates 
the current correlator J E to the gradient of the density 
correlator P E , i.e. a version of Fick's law for the corre- 
lation tensors. This can be achieved by first multiplying 
Eq. i|14f> with q-Q, subsequent integration over the mo- 
menta q and employing the WI, Eq. I|28|) . The resulting 
relation should be evaluated to first order in momentum 
Q at zero external frequency, f2 = (see also the book 
by Case and Zweifel^ for a detailed discussion of this so- 
called Pi-approximation). This procedure leads to Fick's 
law, 



where the renormalization terms A and k are given by 
if 

(27T) 

■d 3 q 



A = J^(q-Q) 2 AG%(0,0)A^ 

3 (g-4)%[ns^(o,o) 



,(o) 



dV 
(2tt) 3 



2tt) 



EG$(0,0W0,0) ■QAG^(0,0)o ( 



(49) 
(50) 



The tensor quantity t 1 (uj) (of dimension inverse length) 

r\u) = is^lSl + AW+ZoV) (51) 

consists of a "standard" term t^ 1 ^) that would also be 
there in electron transport theory, 



rd 3 q 


fd 3 q' 


(2tt) 3 j 


' (2tt) 3 



(«?-Q) 2 AG£(0,0) 



x 7g -, f (0,0)AG|(0,0)a - 1 
d 3 g f d 3 q' 



(2tt) 3 7 (2tt) 3 



(q-Q)tf-Q)AG%(0,0) 



x 7gT (0,0)AG|,(0,0)ar 



(52) 



and two additional contributions g^ 1 1(8>1 and A^' that 
are non-vanishing only in the presence of cither absorp- 
tion or gain, 



A' 3 ' = -iAl?[£± 



(2tt) 3 



(q-Q) 2 nsmo) 



(53) 



+ 



dV 

(2tt) 3 



IIG£ (0,0)7^(0,0) AG^O^r 1 



Similar to the case of the energy transport velocity, any 
discrepancy between the absorption/gain characteristics 
of scatterers and host medium ( "impedance mismatch" 
given by a non- vanishing a^) leads to a renormaliza- 
tion of l~ x . Although the renormalizations A*- 3 ^^) and 
k(lu) originate from isotropic and anisotropic scattering 
effects, respectively, they contribute only in the presence 
of absorption/gain mismatch between scatterer and host 
medium. This is in contrast to earlier results by Livdan 
and Lisyansky 55 for media with frequency independent 
absorbing spheres, due to the fact that these authors in- 
consistently employed the WI without absorption/gain 
when deriving the current-density relation Eq. (|48|l . 

Combining the continuity equation, Eq. (|38|) . with the 
microscopic version of Fick's law, Eq. I|48[) . we are finally 
in a position to solve for the energy density tensor in the 
long-time, long-distance limit 



-g^a + ^))(-^ + D(Q 2 +x- 2 )) 



2 —2 



d 3 qAGm,0) 



(54) 



Q[A + k(lu)] P£ = g^r 1 [1 ® 1 + AM] J%, (48) 



which exhibits a familiar diffusion pole structure. In ad- 
dition, the presence of absorption or gain leads to the 
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appearance of a mass term D(ui)/x~ 2 (lli). This term is 
absent in the work of Livdan and LisyanskyS, although 
they considered absorbing spheres, albeit within a scalar 
model. From Eq. I|54|l the diffusion tensor D(ui) is ex- 
plicitly given by 

D(Lu) = ^{l®t+A(Lu)y 1 £(A + K(Lj)) (55) 

= -v E (uj)e T (uj), 

which establishes the transport mean free path tensor 
£ T (u) as 



3— c P £(A- 



•«H) 



(56) 



Finally, the absorption/gain length tensor x a (u)) explic- 
itly reads 



1 



(57) 



where the absorption mean free path tensor I a is given 
by 



1+A(w))" 



(58) 



Eqs. JSSJl, 0, and represent the central re- 

sults of our investigation and make explicit how the vari- 
ous renormalizations discussed above affect the transport 
quantities of a disordered dielectric medium in the pres- 
ence of linear absorption or gain. In the following section, 
we provide illustrations of these results for a system of 
dilute scatterers. 



VI. NUMERICAL RESULTS 

As alluded to above, any numerical evaluation of Eqs. 
IHHJI7 (G33), an d GUt requires the computation of 
consistent values for the self-energy tensor and the 

irreducible vertex tensor 7~,(Q, 57). For dense systems 
and in the diffusive regime, a Coherent Potential Ap- 
proximation has to be employed for both, the self-energy 
tensor and the irreducible vertex tensor. To date, such 
an approach has been carried out for scalar waves and 
in the presence of point scatterers only&^i. Moreover, 
near the Anderson localization transition, the irreducible 
vertex tensor is expected to exhibit critical behavior and 
more sophisticated schemes such as a self-consistent the- 
ory of localization have to be employed. Again, to date 
such a program has been realized for scaler waves and 
point scatterers only&^Mi. 

However, in realistic strongly scattering electromag- 
netic systems, the scatterers cannot be approximated by 
point scatterers. Instead, they exhibit internal (Mie) res- 
onances which - besides the vectorial nature of the elec- 
tromagnetic radiation itself - greatly complicate the ap- 
proach discussed above. Therefore, we have chosen to 




x = 271 00 / w„ 



FIG. 2: Frequency dependence of the energy transport veloc- 
ity u E units of the vacuum speed of light for a dilute (n — 6% 
by volume) collection of spherical scatterers in air (eh ~ 1) 
with different dielectric constants e s (as indicated). The fre- 
quency unit is luo = 2irc/r where c is the vacuum speed of 
light and r the radius of the scatterers. The pronounced dips 
correspond to the Mie resonances of the individual sphere and 
represent the "dwell time" effect discussed in the text. 



illustrate our results in a simpler system consisting of 
a dilute collection of identical absorbing or amplifying 
spherical scatterers in a host medium without absorp- 
tion or gain. In this case, the independent scatterer ap- 
proximation applies and the self-energy as well as the 
irreducible vertex tensor may be expressed through the 
density n of scatterers and the full t-matrix of such a 
scatterer according tc£2& 



n *&,"?±( w ±) ( 59 ) 

7 W (0; fl ) = (60) 



Explicit expressions for the t-matrix of a spherical scat- 
terers are fairly involved and can be found in a recent 
work by K. Arya et alm&. Consequently, even within this 
approximation, all transport quantities have to be eval- 
uated numerically. 

Note that the integrations with respect to momentum 
of, e.g., certain convolutions of the irreducible Vertex 7, 
which occur in the evaluation of the transport quantities, 
are in general ultraviolet divergent. These divergencies 
are remedied by applying the regularization scheme fa- 
miliar from quantum electrodynamics^, i.e. by subtract- 
ing appropriate mass terms from the integrands such that 
after subtraction, these terms, and sufficiently many of 
their derivatives, do not diverge, thus making the inte- 
grals convergent. Note that the corresponding regular- 
ization masses in general have to be determined numer- 
ically, since the prefactors of the ultraviolet divergencies 
are known numerically only. 

In the following we restrict ourselves to a discussion 
of the numerical results for the energy transport veloc- 
ity tensor. We would like to point out that for isotropic 
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dielectric materials, the energy transport velocity ten- 
sor becomes an isotropic tensor, so that we can further 
restrict our discussion to a single valued, frequency de- 
pendent energy velocity. 

In Fig. |21 we display the frequency dependence of the 
energy transport velocity for a typical concentration of 6 
% by volume of absorption-free scatterers with different 
(real) dielectric constants in an air background. Clearly 
visible are the pronounced dips in the energy transport 
velocity near frequencies that correspond to the single 
scatterer Mie resonances. As the dielectric constant in- 
creases, the Mie scattering becomes stronger and leads to 
correspondingly stronger renormalization of the energy 
transport velocity through the "dwell time" effect men- 
tioned above. These data are, therefore, consistent with 
earlier results for scalar waves^i and vector wavesi^i^. 
This behavior of the energy transport velocity has to be 
contrasted with the corresponding frequency dependence 
of the phase velocity, Eq. (|33|) in Fig. [3J Near the Mie 
resonances, this phase velocity may exhibit values that 
significantly exceed the vacuum speed of light. Since the 
phase velocity in a random medium does not correspond 
to a physical quantity, this behavior is not in conflict with 
the laws of relativity. 

Adding linear gain to the system, dramatically mod- 
ifies the situation. In Fig. 01 we show the frequency 
dependence of the energy transport velocity for the weak- 
est scattering system of Fig. |2 linear gain is added to the 
scatterer. Although the concentration is only n = 6% 
by volume, adding a relatively small negative imaginary 
part to the dielectric constant results in a significant nar- 
rowing of the energy transport velocity resonances. This 
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FIG. 3: Frequency dependence of the phase velocity c P units 
of the vacuum speed of light for a dilute (n = 6% by volume) 
collection of spherical scatterers in air (e^ = 1) with differ- 
ent dielectric constants e s (as indicated). The frequency unit 
is u> = 2nc/r where c is the vacuum speed of light and r 
the radius of the scatterers. Near Mie resonances, this ve- 
locity significantly exceeds the vacuum speed of light. The 
corresponding physical energy transport velocity is displayed 
in Fig. H 
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X = 271 CO / co o 

FIG. 4: Frequency dependence of the energy transport veloc- 
ity u B units of the vacuum speed of light for a dilute (n — 6% 
by volume) collection of spherical scatterers in air (eh ~ 1) 
where the scatterer dielectric constant e 3 exhibits different 
amounts of linear gain (as indicated). The frequency unit is 
uj ~ 2nc/r a where c is the vacuum speed of light and r the 
radius of the scatterers. For increasing gain, the dips in the 
energy transport velocity that are associated with the individ- 
ual scatterers' Mie resonances exhibit considerable narrowing. 
This may be interpreted as a sign for the incipient onset of 
random lasing action in this system. 

gain narrowing may be understood by noting that the 
Mie resonance dips in the transport velocity arise because 
of multiple reflection and interference of light within a 
single scatterer ("dwell time effect"). In the presence of 
gain, the relative importance of long wave paths for the 
interference processes is enhanced, leading to a narrowing 
of the resonance lines, analogous to the narrowing of the 
transmission lines in e Fabry-Perrot filter with increasing 
number of reflections. Increasing the gain to rather unre- 
alistic values shows further narrowing of the resonances 
together with a slight increase in the energy transport 
velocity. It should be noted that the energy transport 
velocity retains physical values for all values of the gain 
that we have considered. 

The opposite behavior occurs when we add linear ab- 
sorption to the weakest scattering spheres of Fig. [21 This 
is illustrated in Fig. \E\ where we display the frequency 
dependence of the energy transport for increasing values 
of linear absorption in this system. In this case, the dips 
in the energy transport velocity that are associated with 
the Mie resonances of the individual scatterers become 
washed out and disappear altogether for sufficiently large 
values of the linear absorption. 



VII. CONCLUSION 

In conclusion, we have presented a microscopic ap- 
proach to the calculation of transport quantities for elec- 
tromagnetic waves propagating in disordered dielectric 



10 



0.8 



0.6 



0.4 




8 s = 2.2 + .01 i 
/ = 2.2 + 1.1 i 
= 2.2 + 5.5 i 
8 S = 2.2 + 12.2 i 



x = 271 00 / 00, 



FIG. 5: Frequency dependence of the energy transport veloc- 
ity v E units of the vacuum speed of light for a dilute (n = 6% 
by volume) collection of spherical scatterers in air (e/j = 1) 
where the scatterer dielectric constant e a exhibits different 
amounts of linear absorption (as indicated). The frequency 
unit is uj — 2irc/r where c is the vacuum speed of light and 
r the radius of the scatterers. For increasing absorption, the 
dips in the energy transport velocity that are associated with 
the individual scatterers' Mie resonances become washed out 
and disappear altogether for sufficiently large absorption, as 
the wave propagation becomes overdamped and interference 
effects become impossible. 



media which exhibit linear absorption or gain. The ef- 
fects of energy conservation and of the violation of en- 
ergy conservation in the presence of absorption or gain 
have been incorporated in the theory by means of an ex- 
act Ward identity. The kinetic equation for light in ran- 
dom media with absorption or gain has been derived and 
solved for the energy density correlation tensor, taking 
fully into account the vectorial nature of electromagnetic 
radiation. In this way we have, to the best of our knowl- 
edge, for the first time derived explicit expressions for the 
energy transport velocity tensor, the diffusion tensor, the 
mean free path tensor, and the absorption length tensor 
in absorbing/emitting media. All these quantities ex- 
perience renormalization due to scattering or impcdcnce 
mismatch between host medium and scatterer material. 
Specifically, we have discussed the case of a dilute sys- 



tem of identical spherical scatterers within the frame- 
work of the independent scatterer approximation. These 
systems exhibit considerable dips in the energy trans- 
port velocity that can directly be traced to a well-known 
"dwell time" effect. However, adding linear gain to the 
scatterers dramatically modifies this situation and signif- 
icant gain narrowing occurs already for relatively modest 
values of the scatterers concentration and the imaginary 
part of their dielectric constant. We interpret this be- 
havior as due to an increase of the relative importance of 
wave paths with long path lengths in the medium, and 
a subsequent, gain- induced enhancement of interference 
effects, analogous to the Fabry-Perrot effect. The oppo- 
site effect occurs when absorption is added to scatterers 
and the resonances in the energy transport velocity are 
washed out and ultimately disappear altogether. In all 
cases, the energy transport velocity retains physical val- 
ues, in contrast to the phase velocity. 

As a systemiatic transport theory, expressing all phys- 
ical quantities ultimately in terms of the single-photon 
sclfenergy and the two-photon irreducible scattering ver- 
tex, our approach may be generalized to include inter- 
ference effects like weak or strong localization, e.g. in 
the spirit of the self-consistent theory of Anderson lo- 
calization that is well-known for electronic systems. To- 
gether with a replacement of the linear gain by a direct 
coupling of transport equations to the scmiclassical laser 
rate equations this may provide a microscopic theory for 
the interplay of optical gain and Anderson localization 
of light. Similarly, the theory may be extended to in- 
clude optically anisotropic materials such as scatterers 
immersed in a liquid crystal or Faraday-active scatter- 
ers. 
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